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ABSTRACT 

Bright submillimetre-selected galaxies have been found to be a rich source of strong 
gravitational lenses. However, strong gravitational lensing of extended sources leads 
inevitably to differential magnification. In this paper I quantify the effect of differ- 
ential magnification on simulated far-infrared and submillimetre surveys of strong 
gravitational lenses, using a foreground population of Navarro- Frenk- White plus de 
Vaucouleurs' density profiles, with a model source resembling the Cosmic Eyelash and 
QSO Jl 148+5251. Some emission line diagnostics are surprisingly unaffected by differ- 
ential magnification effects: for example, the bolometric fractions of [Cn] 158 jam and 
CO (J = 1 — 0), often used to infer densities and ionisation parameters, have typical 
differential magnification effects that are smaller than the measurement errors. How- 
ever, the CO ladder itself is significantly affected. Far-infrared lensed galaxy surveys 
(e.g. at 60 //m) strongly select for high-redshift galaxies with caustics close to active 
galactic nuclei (AGN), boosting the apparent bolometric contribution of AGN. The 
lens configuration of IRAS F10214+4724 is naturally explained in this context. Con- 
versely, submillimetre/millimetre-wave surveys (e.g. 500 — 1400 /im) strongly select for 
caustics close to knots of star formation boosting the latter's bolometric fraction. In 
general, estimates of bolometric fractions from spectral energy distributions of strongly 
lensed infrared galaxies are so unreliable as to be useless, unless a lens mass model is 
available to correct for differential magnification. 

Key words: cosmology: observations - galaxies: evolution - galaxies: formation - 
galaxies: star-burst - infrared: galaxies - submillimetre 



1 INTRODUCTION 

Strong gravitational lensing has been widely exploited for 
studing populations that would otherwise be too faint for 
detailed study (e.g. Smail et al. 1997, Knudsen et al. 2008), 
as well as being one of the very few probes capable of map- 
ping the dark matter distributions in the foreground sys- 
tems (e.g. Gavazzi et al. 2007). Gravitational lenses were 
originally very difficult to discover: for example, the Cosmic 
Lens All-Sky Survey (CLASS) observed 11 685 flat-spectrum 
radio sources, finding 16 gravitational lenses (Myers et al. 
2003). More rapid lens discovery has come from exploiting 
the large public Sloan Digital Sky Survey (SDSS) database. 
The Sloan Lens ACS Survey (SLACS) selected lens candi- 
dates for HST snapshot follow-up from the presence of nebu- 
lar emission lines at a higher redshift than that of the SDSS 
galaxy (e.g. Bolton et al. 2006), resulting in nearly 100 pub- 
lished lenses (Auger et al. 2009). The SLACS lens redshifts 
are z < 0.5 by virtue of their selection criteria, which limits 
the study of the evolving mass profiles of the lenses, while 
the background galaxies are by necessity unobscured blue 
star- forming systems with sufficiently dust-free lines of sight 
through the foreground lens. Many ongoing and future ap- 



plications of gravitational lensing are also limited by sample 
size (e.g. Treu 2010), and so there are strong drivers both 
to finding lenses at higher redshifts and to finding lensed 
systems in greater numbers. 

Several authors predicted that the steep bright-end 
slope of submillimetre source counts, together with the high 
redshifts of submillimetre-selected galaxies (e.g. Chapman 
et al. 2005, Eales et al. 2010) would lead to a strong gravi- 
tational magnification bias (e.g. Blain 1996, Perrotta et al. 
2002, Negrello et al. 2007). Shallow, wide-area submillime- 
tre surveys could therefore be exploited as an extremely ef- 
ficient means of selecting strong lens events, once the con- 
taminant populations of nearby galaxies and blazars are re- 
moved. The millimetre- wave source counts from the South 
Pole Telescope (SPT) were consistent with this prediction 
(Vieira et al. 2010), and the prediction was spectacularly 
confirmed by Negrello et al. (2010), who used an intensive 
and rapid follow-up of the first data from the Herschel Astro- 
physical Terahertz Large Area Survey (if- ATLAS, Eales et 
al. 2010) to find that all five of the first lensed galaxy candi- 
dates are strong gravitational lenses. The selection efficiency 
for strong lenses in bright submillimetre surveys approaches 
an astonishing ~ 100%. The lens redshift distribution is 
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expected to peak at z Z 1, with lensing systems selected 
irrespective of the obscuration in the foreground system. 
The advent of wide-field mapping capabilities at submil- 
limetre and millimetre wavelengths from e.g. Herschel Spec- 
tral and Photometric Imaging Receiver (SPIRE), Submil- 
limetre Common-User Bolometer Array-2 (SCUBA-2) and 
SPT allow the discovery of many hundreds of strong lenses 
(e.g. Vieira et al. 2010, Negrello et al. 2010, Cooray et al. 

2010, Gonzalez-Nuevo et al. 2012). This immediately also 
affords many opportunties for follow-up imaging and spec- 
troscopic diagnostics of the lensed galaxies (e.g. Lupu et 
al. 2010, Cooray et al. 2011, Cox et al. 2011, Frayer et al. 

2011, Gavazzi et al. 2011, Hopwood et al. 2011, Omont et 
al. 2011, Riechers et al. 2011, Scott et al. 2011, Valtchanov 
et al. 2011). 

Gravitational lensing is a purely geometrical effect 
depending on the paths of null geodesies, and is there- 
fore wavelength-independent. Nevertheless, different lines 
of sight have different magnifications in every astrophys- 
ical lens. Therefore if the background source is extended 
with intrinsic colour gradients, the lensed source may have 
colours that are different to the unlensed case. This differen- 
tial magnification could significantly affect broad-band pho- 
tometry and emission line diagnostics. Many strong gravita- 
tional lenses do not have a well- constrained foreground mass 
model, so the level and nature of the differential magnifica- 
tion problem in those sources is unknown. Often, a largely 
unjustified assumption is made, explicity or implicitly, that 
differential magnification effects can be neglected, despite 
the fact that it has long been clear that differential magni- 
fication effects are present and must be corrected for (e.g. 
Blain 1999). 

This paper presents a statistical assessment of the level 
of the differential magnification problem for a wide range 
of photometric and spectroscopic diagnostics of the phys- 
ical conditions in gravitationally lensed infrared-luminous 
galaxies. Section [2] describes the simulation methodology, in- 
cluding the foreground lenses and the adopted background 
source. Sect ion [3] gives the results of the simulations, demon- 
strating where it is and is not appropriate to assume a naive 
demagnification without differential magnification effects. 
Section [H discusses the implications of these results for ongo- 
ing follow-up studies of infrared-selected strong gravitational 
lenses, and the conclusions are summarised in Section[5] This 
paper assumes a concordance cosmology with a Hubble con- 
stant of Ho = 72kms _1 Mpc _1 = 100/ikms" 1 Mpc~\ a 
matter density of Qm = 0.3 and a cosmological constant 
effective energy density of Qa = 0.7. 



2 METHODOLOGY 

2.1 Simulation details 

Differential magnification effects were simulated for a given 
source at a fixed redshift, seen through a constant comoving 
population of lenses. The physical parameters of the lens 
and the source in this case are given in the following subsec- 
tions. The source plane magnification map, summing over all 
images Ei|/Xi(x,2/)|, was calculated using the GRAVLENS soft- 
ware (Keeton 2001). The simulation resolution was 0.01 arc- 
sec at this stage. With the source redshift fixed, the differen- 
tial magnification of the various source components was used 



to calculate observed spectral energy distributions (SEDs). 
A selection at a particular wavelength was then imposed: I 
first assumed that the background source will be considered 
a strong gravitational lens at that observed wavelength when 
the total magnification over all source components exceeds 
a threshold magnification fi > 10. This was then compared 
to the case of more moderate magnifications of 2 < \i < 5 at 
that wavelength. Most of the results presented in this paper 
are for a z = 2 background source, but it was found that 
the results depend only weakly on the background source 
redshift. Differential magnification effects depend on the rel- 
ative angular sizes of the magnification map and the back- 
ground source. As shown in Fig. [I] angular diameter distance 
varies only weakly with redshift at z c± 1 — 3 (see e.g. Serjeant 
2010), while the observed critical radii for singular isother- 
mal spheres depend only very weakly on the source and lens 
redshifts over most of their range (unlike the point mass 
case) . It should therefore not be too surprising that there is 
no strong dependence on the source redshift in these simu- 
lations, notwithstanding changes in rest-frame wavelengths. 
Section I37T1 presents results for a variety of source redshifts. 
Most of the results in this paper are based on lenses selected 
at observed wavelengths of either 60 /mi or 500 /mi, repre- 
senting selection with e.g. IRAS and the Herschel SPIRE 
instrument. The 500 /mi results for a source at z — 2 are 
very similar in practice to simulations at 850 /mi or 1.4 mm, 
representing lens selection with e.g. SCUBA-2 or the SPT; 
this is because all three wavelengths lie well within the red- 
shifted Ray leigh- Jeans tail. Section T3.7I discusses the depen- 
dences on observed wavelength in more detail. 

A higher resolution (0.001 arcsec) simulation was then 
run for the specific case of a lens geometry based on that 
of IRASFSC 10214+4724, with a single lens at a redshift of 
z — 0.9 and a source at z = 2.286. Finally, calculations were 
made of flux- or luminosity-limited samples in the context 
of a particular source count model. 



2.2 The foreground lenses 

The foreground lensing galaxy or galaxies were simulated 
following SLACS (Bolton et al. 2006). The assumed density 
profile was the sum of a Navarro- Frenk- White (NFW) pro- 
file and a de Vaucouleurs' profile, which has an azimuthally- 
averaged surface density similar to a singular isothermal el- 
lipsoid (e.g. Gavazzi et al. 2007). Spiral galaxy discs present 
only 10 — 20% of the lensing cross-section of elliptical galax- 
ies (e.g. Keeton & Kochanek 1998, Moller et al. 2007) so 
spiral discs are neglected in this analysis. Also neglected 
for simplicity is the lensing by galaxy groups and clus- 
ters, in which it may be rarer to find subarcsecond-scale 
strong magnification gradients intercepting the background 
sources. Nonetheless, the existence of giant optical arcs im- 
plies that galaxy clusters can sometimes induce strong mag- 
nifications (e.g. Bayliss et al. 2011, Wuyts et al. 2012). In 
these situations, one might expect differential magnification 
to be relevant. 

The NFW profiles of the elliptical galaxy haloes in this 
paper had a virial mass equal to the mean SLACS lens mass, 
Mvir = 14/i -1 x 10 12 M©, and a concentration parameter c 
given by 
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Figure 1. The Einstein radius in arcseconds against lens redshift 
for a singular isothermal sphere with velocity dispersion equal 
to the SLACS mean of 230 km s -1 (upper curves, Bolton et al. 
2008), and for comparision a one solar mass point lens (lower 
curves), at source redshifts of 0.5, 1, 2, 4 and 8 (dotted, dash- 
dot, dashed, long dashed and full lines respectively). Note that 
the Einstein radius for the singular isothermal sphere has a much 
weaker dependence on the lens redshift than the point mass case, 
until the lens approaches the source. As a further comparison, 
the dotted line shows the observed angular size of a fiducial one 
parsec proper length. Note that over z ~ 0.5 — 3, the observed 
angular sizes are relatively weakly dependent on redshift. The 
weak dependence on redshift of both the source sizes and lens 
critical radii imply the results in this paper will be only weakly 
dependent on the background source redshift, notwithstanding 
differences in rest-frame wavelengths. 
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(Hoekstra et al. 2005, correcting a typographical error in the 
h dependence, e.g. Bullock et al. 2001, Gavazzi et al. 2007). 
The virial radius was set to R v [ r = 466//ikpc, consistent 
with 
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where p c is the critical density and A = 119. For the de 
Vaucouleurs' profile, the effective radius was set to the mean 
of the SLACS lenses (Gavazzi et al. 2007), 5.3/i _1 kpc. 




[Cji] Cirrus + 
cool CO 



Figure 2. A sketch, not to scale, of the emission regions in the 
model background source described in Section [2?3l The GMCs are 
embedded in larger warm CO regions. 



For the case of a single lens, the configuration 
was modelled on the lensed hyperluminous galaxy IRAS 
F10214+4724 (Rowan- Robinson et al. 1991, Broadhurst & 
Lehar 1995, Graham & Liu 1995, Serjeant et al. 1995, Eisen- 
hardt et al. 1996), with a background source at a redshift of 
z — 2.286 (e.g. Brown & Vanden Bout 1991) and a lens at 
z = 0.9 (e.g. Serjeant et al. 1995, 1998). The ellipticity of the 
lens, denned as one minus the axis ratio, was set to e = 0.21, 
equal to the median value for SLACS lenses (Bolton et al. 
2008). No external shear was included as only one angle de- 
pendence is necessary to reproduced the SLACS lens mor- 
phologies. Under these assumptions, the critical (Einstein) 
radius is 1.18 arcsec, close to the median value of 1.17 for 
the SLACS singular isothermal ellipsoid models (table 5 of 
Bolton et al. 2008). 

For the case of a population of lenses, a constant comov- 
ing density of these SLACS-like lenses was assumed, whose 
mass profiles do not evolve except in the variation of concen- 
tration parameter (equation [T]) . No assumptions were made 
about the comoving density of background sources; instead, 
the differential lensing statistics were found as a function of 
source redshift. 



2.3 The background source 

2.3.1 SMG model 

I modelled a submillimetre-selected galaxy (SMG) -like ex- 
tended background source with a structure resembling that 
of the Cosmic Eyelash (Swinbank et al. 2010) and the 
z = 6.42 quasar QSO J1148+5251 (e.g. Walter et al. 2009). 
The objective is to calculate the differential magnification ef- 
fects on common emission-line diagnostics (CO, [Cn], H2O, 
etc.) and bolometric fractions (AGN, star formation, cir- 
rus) so the model background galaxy has spatially localised 
regions and associated SEDs for each of these emission com- 
ponents. Fig. [2] shows a sketch (not to scale) of the compo- 
nents of this simulated galaxy. The source has the following 
features: 

(i) An AGN with a bolometric fraction of 0.3, with an 
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SED from Nenkova et al. (2008a, b), modelled as a circular 
region with 0.1 kpc radius. The results are not strongly de- 
pendent on the AGN torus model chosen, as long as its SED 
peaks in the rest-frame mid-infrared. The selected model has 
an AGN and dust torus, incident AGN spectum AA01, an 
angular distribution width a = 30°, a Gaussian angular dust 
distribution with an equivalent width of 45°, a radial torus 
thickness of 100 sublimation radii, an average of 14 clumps 
along a radial equatorial ray, a cloud distribution with a 
radial dependence of r -20 , and a visual optical depth of 
individual clumps of tv = 100. The torus is simulated as 
viewed from an inclination angle of 50° . 

(ii) An extended cirrus region with a 2.5 kpc major axis, 
an ellipticity of e = 0.4, a bolometric fraction of 0.2, and 
a cirrus SED from Efstathiou, Rowan- Robinson & Sieben- 
morgen (2000). 

(iii) Four extreme giant molecular cloud (GMC) systems, 
as seen in the Cosmic Eyelash, with a total bolometric frac- 
tion of 0.5. The GMCs (assumed identical) are placed ran- 
domly into the cirrus region, and have a starburst SED taken 
from the models of Efstathiou, Rowan- Robinson & Sieben- 
morgen (2000), with an age of 1.7 Myr and an optical depth 
of rv = 100. They are modelled as 50 pc radius circular 
clouds. 

(iv) A circular old stellar population with a 2 kpc half- 
light radius and a de Vaucouleurs' profile. The 2 Gyr-old el- 
liptical SED from the Spitzer Wide- Area Extragalactic Sur- 
vey (SWIRE) library (Polletta et al. 2007) is adopted for 
this component, though for this study it is only necessary 
that some relatively AGN-free observational proxy is avail- 
able (e.g. Spitzer 3.6 /mi fluxes) and that it is possible to 
subtract the foreground lens (e.g. Hopwood et al. 2011). 

(v) A circular [Cn] emitting region offset 0.6 kpc from the 
AGN with radius 0.75 kpc, following the precedent of SDSS 
QSO Jl 148+5251. Similarly, the lensed star- forming galaxy 
MIPS J142824.0+352619 has been argued to contain a [Cn]- 
emitting region extended over kpc scales (Hailey-Dunsheath 
et al. 2010). 

(vi) HCN and HNC emission regions are modelled as two 
point sources separated by 80 pc (see Section T3.3p . 

(vii) Two H2O emission regions: a circular "warm" region 
with a radius of 120 pc, and an elliptical "cool" region with 
a major axis of 1 kpc and an ellipticity and position angle 
matching the cirrus component. This follows the precedent 
of Mrk231 (van der Werf et al. 2010, Gonzalez- Alfonso et 
al. 2010). 

(viii) A cool molecular gas component tracing the cir- 
rus. The CO emission spectrum was calculated using the 
RADEX software (van der Tak et al. 2007) using the Large 
Velocity Gradient (LVG) approximation, assuming it is com- 
prised of clumps in a background temperature of 8.97 K (i.e. 
2.73(1 + z) K at z — 2.286), with a kinetic temperature of 
20 K, an H2 density of 100 cm -3 and a CO column density 
of 10 14 cm -2 for each cirrus cloud. 

(ix) A warm molecular gas knot with 0.4 kpc radius cen- 
tred on each GMC. Following the precedent of SDSS QSO 
Jl 148+5251, these warm components in total contribute 1/3 
of the total CO (J = 3 — 2) flux. The CO emission of these 
knots was calculated using the RADEX (van der Tak et al. 
2007) LVG model assuming they are comprised of clumps 
in a background temperature of 8.97 K (i.e. 2.73(1 + z) K at 
z — 2.286), with a kinetic temperature of 45 K an H2 density 



of 10 4 cm -3 and a CO column density of 10 18 cm -2 . There 
is already good evidence that the higher J transitions of CO 
in submillimetre-selected galaxies are much more spatially 
localised than CO(J =1-0) (e.g. Ivison et al. 2011). 

This model source will be referred to as the SMG model. 
In the simulations with 0.01 arcsec resolution, the max- 
imum monochromatic magnification of a 500 /mi-selected 
lens at z — 2 for a constant comoving density of lenses 
is ji = 33. Restricting this to configurations with magni- 
fications \i > 10, the median (mean) magnification is 11.5 
(12.0). At an observed wavelength of 60 /mi the flux is dom- 
inated by much smaller scale features such as the AGN, so 
the corresponding numbers are higher: the median (mean) 
magnification of fi > 10 sources is 18.5 (21.9), while the 
maximum monochromatic magnification is 840. The config- 
urations that approach this (possibly unphysical) maximum 
are nonetheless extremely rare, and the results in this pa- 
per are unaffected by excluding all configurations with mag- 
nifications of e.g. \i > 100. The higher angular resolution 
simulation discussed in Section 13.61 has a more moderate 
maximum magnification. 



2.3.2 Filling factors 

The filling factor of emission line gas in SMGs is far higher 
than in local galaxies. For example, Carilli et al. 2010 derived 
a filling factor of 0.5 in the low excitation gas component of 
their z — 4.05 SMGs, and 0.13 in its higher excitation gas. 
These values are an order of magnitude or more larger than 
seen in local galaxies (e.g. McKee & Ostriker 2007). In effect, 
the clumping of the emission line gas in high-redshift star- 
forming galaxies creates a sparse sampling of the morpho- 
logical distributions defined above. Since the surface filling 
factor is high, there are unlikely to be problems caused by 
small number statistics of clumps. The objective is to simu- 
late a high-redshift galaxy, and therefore in the SMG model 
the emission line gas is treated as spatially continuous. Any 
further clumping beyond that already specified is neglected. 



3 RESULTS 

3.1 Emission line diagnostics: [Cn] and [Oi] 

Fig. [3] shows a simulated emission line diagnostic diagram 
for a 500 /mi-selected SMG model galaxy, lensed with a mag- 
nification of \i > 10 and of 2 < \i < 5, with [Cii]/Lfir plot- 
ted against CO (J = 1 — 0)/Lfir. This diagram has been 
used to estimate the density and ionizing radiation envi- 
ronment in star-forming galaxies (e.g. Stacey et al. 2010, 
Valtchanov et al. 2011). The 500 /mi luminosity is taken as 
a proxy for the far-infrared (FIR) luminosity Lfir in the 
differentially-magnified simulated galaxies. It is clear that 
differential magnification has surprisingly little effect on this 
diagnostic plot for lensed galaxies. For nearly all lens config- 
urations, differential magnification confines the source to a 
line almost parallel to the constant density contours. How- 
ever, the vertical positioning of this track was found to de- 
pend on the random positioning of the GMCs in the source 
model. 

If one chooses to use the CO (J = 4 — 3) line as a proxy 
for the CO (J =1 — 0) line, the smaller emitting regions of 
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Figure 3. Simulated emission line diagnostic diagram for a 500 /xm-selected gravitational lens with a source at z = 2; very similar results 
were found for sources at other redshifts. The underlying value of the line/FIR ratios is taken to be identical to the H- ATLAS galaxy 
ID81 (Valtchanov et al. 2011), while the differentially magnified simulated data points are shown as red contours near the ID81 data 
point (contours at 68%, 95%, 99% and 99.9%). The observed data are taken from the compilation in Stacey et al. 2010 and Valtchanov 
et al. 2011. Also shown are densities and interstellar far-ultraviolet flux predictions from Stacey et al. 2010; densities are in units of cm -3 
and the ionising flux Fpuv is parameterised as Go = Fpuv/(l-6 x 10 -6 Wm -2 ). The left-hand panels show the use of the CO(J =1 — 0) 
line, while the right-hand panels show the use of CO(J = 4 — 3) as a proxy for CO(J = 1 — 0). The lower panels show results for moderate 
magnifications of 2 < \i < 5, while upper panels show high magnification systems of \i > 10. In all four cases, nearly all the contours 
follow a track that is almost parallel to a constant density diagnostic curve. The vertical positioning of this track was found to depend 
on the source structure. 



CO (J = 4 — 3) slightly changes the differential magnification 
effects, as shown also in Fig. [3] Table [T] gives the dispersions 
in the logarithmic bolometric fractions of the emission lines 
in Fig. [3j The dispersion in log 10 CO/Lfir is smaller than 
the measurement errors of if- ATLAS galaxy ID81, while the 
standard deviation of log 10 [Cn]/Z/FiR is comparable to the 
measurement errors for ID81 (Valtchanov et al. 2011). 

The SHINING survey (Survey with Herschel of the 
ISM in Nearby Infrared Galaxies, e.g. Sturm et al. 2010) 
has found that the [Ol] 63 /xm emission region in M82 is 
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g 10 (CO(J = 1 - 0)/FIR) 


0.16 


g 10 (CO(J = 4-3)/FIR) 


0.15 


gio([Cn]/FIR) 


0.17 



0.10 
0.10 
0.19 

Table 1. Dispersions in the logarithmic bolometric fractions of 
key emission line diagnostics from Fig. [3] for the stated magnifi- 
cations ji. 
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Figure 4. Simulated CO Spectral Line Energy Distribution (SLED) for a z = 2 gravitational lens, selected at an observed wavelength 
of 60 /im (left) and 500 /im (right), using the SMG model as the source. Very similar results were obtained for other source redshifts. 
The upper panels show the SLEDs for high-magnification lenses with \i > 10, while the lower panels show the SLEDs for more moderate 
magnifications of 2 < \i < 5; only a random sparse sample of 2 (10) percent of the lens configurations are shown in the moderate (strong) 
magnification case, for clarity. The unlensed cool component is shown as the blue line, while the unlensed warm component is a red line. 
The magenta line and large symbols show the sum of the unlensed compoents. The total lensed CO SLEDs are shown as dotted lines, 
with the dots spaced randomly along each line to avoid overlapping points. In each case, one particularly favoured configuration leads to 
multiple superimposed dotted lines (i.e. the black line). All SLEDs have been normalised to the CO( J = 1 — 0) point, to remove the the 
overall effect of magnification and isolate the effects of differential magnification. 



broadly co-spatial with that of the [Cn] 158 /im emission 
on 100 — 200 pc scales (Contursi et al. 2010, Sturm et al. in 
preparation) . If this is also true of high-redshift star- forming 
galaxies, the density-sensitive [Oi]/[Cii] line ratio should be 
insensitive to differential magnification. 



3.2 Emission line diagnostics: CO 

Cox et al. (2011) find the z = 4.24 lensed SMG ID141 has 
a CO Spectral Line Energy Distributions (SLED) similar in 
shape to that of the z = 4.05 SMG GN20 (Carilli et al. 



2010) and in the z = 2.5 lensed galaxy SMM J16359+6612 
(WeiB et al. 2005), though cooler than that of the z = 5.3 
millimetre-galaxy AZTEC-3 (Riechers et al. 2010). Could 
differential magnification skew the CO SLEDs in the lensed 
galaxies? 

Fig. [3] already suggests the CO ladder is distorted due to 
the differing spatial extents of the CO emission line regions 
in our source model. In this case, the effect on the physical 
interpretation is much stronger, because of a much smaller 
dynamic range used for the diagnosis of physical conditions. 

Fig. g] shows the SLEDs for the lensed SMG model 
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Figure 6. Apparent bolometric fractions for the simulated z = 2 SMG model source in Section I2.31 gravitationally lensed by the 
lens population described in Section T2.2I then selected at 60 fim (left) or 500 jum (right). The upper panels show the results for high 
magnification lenses, /i > 10; the lower panels show moderate magnifications lenses, 2 < /i < 5. The large red dot represents the 
underlying unlensed source, while the small dots represent the differentially magnified source. The data are represented as a ternary 
diagram, i.e. the perpendicular distance from any side to the opposite corner is proportional to the bolometric fraction, with 0% on the 
side and 100% at the corner. Each corner is labelled with the 100% bolometric fraction at that point. Only a random sparse sample of 2 
(10) percent of the configurations are plotted in the lower (upper) panels, for clarity. 



source. The SLEDs have been normalised to the CO (J = 
1 — 0) point to remove the overall magnification effects and 
isolate the differential magnification effects. Clearly the rel- 
ative proportions of the cool and warm molecular gas com- 
ponents vary strongly, according to the random differential 
magnification, even for moderate magnification lenses. Both 
the 500 /xm-selected lenses and 60 /mi-selected lenses suffer 
strong differential magnification effects, though the effects 
are not identical. There is also a higher-probability SLED 
(the darker bands in Fig. 3} within the overall dispersion 
that was found to depend on the specific random configu- 
ration of the GMCs. The warm molecular gas component 
is in close physical proximity to the star-forming knots in 
this model, but this does not necessarily guarantee that the 
emission from the warm gas component is boosted relative 



to that from the cool component, even when selecting at an 
observed wavelength of 500 /mi. The la dispersion in the ap- 
parent CO (J = 6 — 5) line (once normalised to CO (J =1 — 0) 
as in Fig. 3J is 32% of the mean value for 60 /xm-selected 
strong lenses (/i > 10), or 20% of the mean value for 500 /xm- 
selected strong lenses. Even at more moderate magnifica- 
tions of 2 < \i < 5, the corresponding dispersions are 29% 
at 60 /im and 30% at 500 /mi. 



3.3 Emission line diagnostics: HNC, HCN, H 2 

Omont et al. (2011) found the bolometric fraction of an H2O 
line in the z = 2.3 lensed SMG SDP 17b to be higher than 
that in Mrk231, suggestive of a luminous AGN. Van der 
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Figure 5. Simulated differential magnification of the H2O emis- 
sion regions in the z = 2 SMG model discussed in the text, for a 
500 /im-selected gravitational lens with a magnification /i > 10. 



Werf et al. (2011) detected four rotational H2O transitions 
in the z = 3.9 lensed quasar APM 08279+5255, deriving a 
gas temperature of 105 =b 21 K. Lupu et al. (2011) reported 
early results from high density gas tracers in lensed SMGs, 
such as HCN, HNC HCO + and 13 CO. Are these observations 
sensitive to the effects of differential magnification? 

The J — 1 — HNC and HCN transitions at rest-frame 
wavelengths of around 3 mm have very similar critical densi- 
ties (3 x 10 6 cm -3 ) and so there is a reasonable expectation 
for these lines to be co-spatial. The HCNiHNC flux ratio is 
a good diagnostic of whether the emission line gas is from 
a photodissociation region (PDR) or an X-ray dissociation 
region (XDR), e.g. Loenen et al. 2007. The HCN:HCO + and 
HNC:HCO + ratios, conversely, are sensitive to densities in 
PDRs, since HCO + has a slightly lower critical density of 
3x 10 5 cm" 3 . Aalto et al. (2009) find that the J = 3-2 HCN 
and HNC emission in Arp 220 are both confined to a region 
< 350 pc wide, with most within the central 30 pc. The HCN 
and HNC emission has therefore been pessimistically mod- 
elled in the SMG model source as point sources separated by 
80 pc. The simulations confirmed that the differential mag- 
nification effect are essentially negligible between these two 
emission regions. 

A similar optimistic situation exists for moderate mag- 
nification in H2O lines. Following the precedent of Mrk231, 
the H2O emission is modelled as two components, one with 
a radius of 120 pc, and the other with a major axis of 
1 kpc (Section 12. 3p . These components contribute compa- 
rable flux in Mrk231 at upper energy levels below 200 K, 
while the warm component dominates at higher tempera- 
tures. For 500 /im-selected lenses with moderate magnifica- 
tions of 2 < \i < 5, the simulations show almost no differen- 
tial magnification between the two components. Very similar 
results were obtained for other source redshifts. 

However, the situation is not as optimistic for H2O lines 
at magnifications \i > 10. These configurations have strong 
differential magnification effects between these two compo- 
nents, shown in Fig. [5] This figure also shows a tail to high 
differential magnifications. A similar tail is also present in 
the magnification ratio of warm H2O to 500 /mi emission, 
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Figure 7. Histogram of apparent AGN bolometric luminosity 
fractions in the simulated z = 2 SMG model source from Sec- 
tion 12.31 and a lens population as described in Section 12.21 for 
lenses with magnification \i > 10 selected at an observed frame 
of 60 /im and 500 /im. The underlying AGN bolometric fraction is 
0.3 in this model (Section 12. 3|) . On average, the 60 /im selection 
over-represents configurations in which the AGN is close to a caus- 
tic and has a consequent strong differential magnification. Con- 
versely, 500 /xm selection more often has GMCs with the greater 
magnification. Nevertheless, note that many lens configurations 
that satisfy the 60 /im lens selection would also satisfy the 500 /im 
lens selection in this simulation. 



and in that of the warm H2O to GMC emission. The tail is 
an alternative interpretation of the H20(2o2 — In) detection 
in the lensed H erschel- ATLAS galaxy SDP.17b (Omont et 
al. 2011), since in Mrk231 about two-third of the luminosity 
in this line comes from this compact warm component. 

The ratios of H2O to high- J CO transitions have also 
been argued as a useful diagnostic of PDR versus XDR con- 
ditions (e.g. Gonzalez- Alfonso et al. 2010, van der Werf et al. 
2010). In highly magnified lensed systems this will clearly be 
contingent on the degree that the emissions can be treated 
as co-spatial. 



3.4 Apparent bolometric fractions 

Broad-band SEDs have long been used to derive the relative 
bolometric fractions in galaxies of AGN and star formation, 
via radiative transver models (e.g. Rowan- Robinson 1980, 
2000; Rowan- Robinson & Crawford 1989; Granato, Danese 
& Franceschini 1996; Green & Rowan- Robinson 1996; Silva 
et al. 1998, 2011; Efstathiou, Rowan- Robinson & Siebenmor- 
gen 2000, 2006; Popescu et al. 2000; Misiriotis et al. 2001; 
Farrah et al. 2002; Verma et al. 2002; Tuffs et al. 2004; 
Mollenhoff, Popescu & Tuffs 2006; Popescu et al. 2011). 
However, it is rare that differential magnification is explicitly 
incorporated into SED fits of lensed star-forming galaxies 
(e.g. Deane et al. 2011). 

It is very clear from the lensed SMG model that differ- 
ential magnification has a strong effect on the apparent bolo- 
metric fractions of the continuum components in the model 
source. Fig. [6] shows the effect of differential magnification 
on these components in the SMG model source, expressed 



Strong biases in infrared- selected gravitational lenses 9 







0.0 0.5 1.0 1.5 2.0 2.5 

Mstars/MSFR 

Figure 8. The shaded histogram (bin size 0.05) shows the rel- 
ative magnification factor of the old stars component with that 
of the star forming regions in the z = 2 SMG model source with 
magnifications a > 10. The unshaded histogram (bin size 0.02) 
gives the result when using the 500 /im luminosity as a proxy 
for the star formation rate. Both histograms are for a simulated 
500 /im-selected lens. 

as a ternary diagram. At each apex, the bolometric fraction 
of that component is 100%. Along each side, the bolometric 
fraction of the opposite component is 0%. For example, in 
the bottom-left corner, the bolometric fraction of the cirrus 
component would be 100%, while all along the right-hand 
side the cirrus bolometric fraction is 0%. The bolometric 
fraction is proportional to the perpendicular distance from 
a side. 

The large red dots in Fig. [6] represents the underlying 
bolometric fractions in the unlensed source, from Section l231 
The small dots represent the simulated magnified source, for 
60 /mi-selected lenses and for 500 /mi-selected lenses. The 
cirrus does not have a strong differential magnification, be- 
cause its physical size is large compared to the caustic sizes 
in the source plane. However, the GMCs and the AGN are 
small, and therefore subject to strong differential magnifica- 
tion effects. 

In the case of high magnifications (/j, > 10), these dif- 
ferential magnification effects are very large. It is clear from 
Fig. [6] that without a lens model to correct for differential 
magnification in strongly lensed systems, the observed bolo- 
metric fractions of AGN, starburst and cirrus are so unreli- 
able as to be useless for diagnosing the underlying bolomet- 
ric fractions. 

A curious feature of Fig. [6] are the unoccupied regions 
for strong lenses (/i > 10). The spatial extent of the cirrus 
limits its maximum magnification, which corresponds to a 
diagonal threshold parallel to the right-hand side. However, 
there are clearly other thresholds also at work: for example, 
there are no simulated lens configurations with a maximum 
cirrus bolometric fraction and an AGN bolometric fraction 
of ~ 0%. This reflects the specific spatial configuration of the 
background source, i.e. there are no configurations in which 
a GMC is close to a caustic while the AGN is simultaneously 
distant from all caustic curves. 

The different distributions in Figs. [6] demonstrate that 



the selection wavelength imposes a strong bias on the ob- 
served apparent bolometric fractions at high magnification 
factors. These differences between the observed- frame se- 
lection wavelengths are clearer in Fig. [71 which shows the 
AGN bolometric fraction for the two selection wavelengths 
for /x > 10 lenses. Selecting at 60 /mi is sampling < 20 /mi 
in the rest-frame, so one may expect 60 /mi selection to 
favour configurations the AGN torus is close to a caustic 
in the source plane. This is indeed what is seen in Fig. [7] 
This also naturally explains the lens configuration of IRAS 
F10214+4724, in which a cusp appears to lie across the AGN 
dust torus (Lacy, Rawlings and Serjeant 1998). A similar 
suggestion was made for other lensed ultraluminous galax- 
ies by Blain (1999). Conversely, 500 /mi selection favours the 
differential magnification of the most luminous small sources 
at 500 /mi, which are the GMCs. Note however that in this 
simulation, many lens-source configurations that would sat- 
isfy the 60/mi-based lens selection, would also satisfy the 
500/mi-based lens selection; it is simply that there are ad- 
ditional configurations that are are more common in the 
500 /mi case (Fig. 0). 

For moderate magnification lenses (2 < \i < 5), the ob- 
served AGN fraction is an unbiased estimator of the under- 
lying value. In these simulations, the apparent AGN bolo- 
metric fraction for both 60 /mi-selected and 500 /mi-selected 
galaxies is roughly normally distributed about the underly- 
ing value of 0.3, with a dispersion of 0.03. 



3.5 The galaxy "main sequence" 

The suggestion has been made that star-forming galaxies 
have a characteristic specific star formation rate, which 
evolves from z = — 2 and then plateaus (e.g. Noeske et 
al. 2007, Elbaz et al. 2007, Daddi et al. 2007, Stark et al. 
2009, Gonzalez et al. 2010). The scatter in this relationship 
is around a factor of two or more. There has been some 
discussion as to whether submillimetre-selected galaxies are 
typical of the population as a whole, or whether they are 
outliers in this relation (e.g. Tacconi et al. 2006, Hainline et 
al. 2011, Michalowski et al. 2012). Could differential magni- 
fication drive a lensed submillimetre galaxy into appearing 
as an outlier? 

This effect would be strongest at high magnifications, 
so the \i > 10 case will be discussed first. Selecting the SMG 
model lenses at 500 /mi, and assuming that the old stellar 
population can be observed at a wavelength at which AGN 
contributions are negligible and that the foreground lens can 
be subtracted (e.g. Hopwood et al. 2011), the observed flux 
ratio of the old stellar population to that of the star- forming 
regions varies by ±0.26 dex (lcr), as shown in Fig. [8] (This 
also assumes that one can successfully address problems of 
the choice of initial mass function, evolutionary synthesis 
models, and star formation histories in estimating the stellar 
masses; e.g. Michalowski et al. 2012.) The median differen- 
tial magnification ratio is 0.86. Using the 500 /mi emission 
itself as a proxy for the star formation rate, the dispersion 
reduces to ±0.07 with a median differential magnification 
ratio of 0.77 (Fig.®. 

For more modest magnifications (2 < ji < 5), the ef- 
fects are much smaller. The distributions are centred on 
Mstars/MSFR = 1 with a dispersion of 0.14, or 0.06 if using 
500 /im luminosity as a proxy for the star formation rate. 
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Figure 9. A selection of results for simulated high magnification (fi > 10) lensing configurations for the source and lens redshifts of 
IRAS FSC 10214+4724. Symbols as in Figs. H [6] and 



There are two lessons to be drawn from this. Firstly, 
the specific star formation rates can be misleading if the 
star formation rate proxy is also the flux at the selection 
wavelength. This is closely related to the fact that bolomet- 
ric fractions from the SED shape can be misleading, as noted 
in the previous section. Secondly, the dispersion induced by 
differential magnification in this case is comparable to the 
rather broad dispersion in the galaxy "main sequence". A 
large offset from this correlation would therefore be unlikely 
to be attributable to differential magnification. 



3.6 Results for an IRAS FSC 10214+4724-type 
lens 

How sensitive are the simulations to the angular resolution? 
The simulation was repeated a resolution of 0.001 arcsec, i.e. 
ten times finer than the lens population above, for a lens at 
a fixed source redshift of z = 2.286 and a lens at z — 0.9, 
modelled on the example of IRAS FSC 10214+4724. The ori- 
entation of the lens ellipticity was fixed at an arbitrary value 
of 7r/4 relative to that of the source. The maximum magnifi- 
cation at an observed wavelength of 500 /xm was 17, while at 



an observed wavelength of 60 /mi it was 52. Fig. [9] shows a 
selection of results for this case. The results are qualitatively 
similar to the z — 2 lens case with a constant comoving pop- 
ulation of lenses; the peaks in the AGN bolometric fraction 
histograms are due to particular alignments of source struc- 
ture with cusps in this model. No significant quantitative 
differences in Fig. [9] were found with an otherwise identical 
simulation at 0.01 arcsec resolution, justifying the resolution 
choice of the lens population simulations above, despite the 
fact that the smallest source structures had been treated as 
de facto point sources. More results from this lens configu- 
ration are found in Serjeant (2011). 



3.7 Results for other wavelengths and source 
redshifts: ALMA and SMA 

How does differential magnification depend on source 
redshift and observed wavelength, and how can 
submillimetre/millimetre-wave continuum photometry 
and interferometry illuminate the bolometric fractions? 

The discussion up to this point has been restricted to a 
source redshift of z — 2.0, with the expectation from Fig.[T] 
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that there should not be a strong dependence on the source 
redshift, with the exception of a possible dependence on the 
rest-frame wavelengths being sampled. Confirmation of this 
can be found in Fig. [10] in which the apparent AGN bolo- 
metric fraction is plotted for source redshifts of z = 1 and 
z = 4. 

Most of the results in this paper for 500 /mi-selected 
galaxies apply equally to 1.4mm-selected galaxies, since 
both wavelengths usually lie well within the Ray leigh- Jeans 
tail. However, at sufficiently high redshifts, this situation 
must change. Fig.[TT]presents a comparison of the monochro- 
matic magnifications at 500 /xm and 1.4 mm for source red- 
shifts of z = 1, z = 2 and z = 4. Only at the high- 
est redshifts and high magnifications are there apprecia- 
ble differences in the monochromatic magnification fac- 
tors. A corollary is that e.g. 870 /xm Sub-Millimeter Ar- 
ray (SMA) follow-ups of 500 /xm-selected lenses will be ac- 
curately recovering the 500 /mi magnifications and mor- 
phologies, provided the 500 /xm passband is not close to 
the peak of the SED, or that the magnification factor is 
not large. Similarly, SMA mapping will almost always re- 
cover the 1.4 mm morphologies and magnifications. ALMA 
mapping of submillimetre/millimetre- wave-selected strongly 
lensed galaxies will also recover the morphologies and mag- 
nifications at the selected wavelength, provided that both 
the ALMA band and the selection wavelength are on the 
Ray leigh- Jeans tail (or that the ALMA band and the selec- 
tion wavelengths are similar). 

However, it does not necessarily follow that single-band 
photometry is a good tracer of any particular component 
of the background galaxy. Fig. [12] correlates the monochro- 
matic magnification factors at selected wavelengths with the 
magnifications of the AGN, starburst and cirrus compo- 
nents. As one might expect, mid-infrared rest-frame pho- 
tometry is a good tracer of the magnification of the AGN 
component. The magnification of the cirrus component is 
well-traced by photometry on the Ray leigh- Jeans tail. The 
starburst component, however, appears not to be well-traced 
by any single monochromatic photometric measurement. 
In order to decouple the starburst and cirrus components, 
one must either obtain more extensive coverage of the FIR 
to submillimetre/millimetre- wave SED, or obtain resolved 
imaging to decouple the FIR-luminous components. 



3.8 Source count model predictions 

One potential bias not accounted for in the discussion so far 
is the effect of the slope of the luminosity function. As the 
slope steepens, the magnification bias is stronger. Might one 
expect more pathological lens configurations as one samples 
populations progressively brighter than the break luminosity 
L*? 

Lapi et al. (2011) recently presented constraints on 
the evolving 100 /xm rest-frame luminosity function of 
submillimetre-selected galaxies. Their z — 2.4 — 4 luminosity 
function, and their toy model that fits it, are both consistent 
with a Schechter function with a faint-end slope of —1.46. 
This luminosity function was assumed to hold at z = 4. 
Submillimetre galaxies were sampled from this distribution 
and magnified by the constant comoving density of lenses 
described above. The results are shown in Fig. [13] It is clear 
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Figure 10. The apparent AGN bolometric fractions for a source 
redshift of z = 1.0 and z = 4.0, for /i > 10 lensed galaxies selected 
at 60 /im (grey) and 500 /im (open). Note the similarity to Fig.[7] 



that the qualitative results of the magnification simulations 
in Section \3. 4 1 are unaffected. 



4 DISCUSSION 

The obvious underlying principle of these simulations is that 
spatially-compact regions are prone to differential magnifi- 
cation, particularly in the more strongly- magnified systems, 
while spatially-extended regions are less prone to this bias. 
For the specific SMG model galaxy and foreground lens pop- 
ulation discussed in this paper, this underlying principle 
leads to several specific conclusions with clear implications 
for the follow-ups of infrared-selected strong gravitational 
lens systems. 

It is particularly clear that for high-magnification lenses 
(/i > 10), the observed selection wavelength has a strong 
effect on the distribution of relative bolometric fractions. 
Radiative transfer models are an extremely effective and 
well-established route to decomposing the relative bolomet- 
ric fractions of ultraluminous galaxies (e.g. Rowan-Robinson 
1980; Rowan- Robinson & Crawford 1989; Green & Rowan- 
Robinson 1996; Granato, Danese & Franceschini 1996; Silva 
et al. 1998, 2011; Efstathiou et al. 2000; Popescu et al. 2000, 
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Figure 11. The monochromatic magnification factors at ob- 
served wavelengths of 350 /im, 500 /i-m, 850 /im and 1.4 mm, for 
source redshifts of z = 1, z = 2 and 2 = 4 and lens redshift of 
half the source redshift. Note that the magnifications are essen- 
tially identical until the passband moves off the Rayleigh- Jeans 
tail. 



2011; Misiriotis et al. 2001; Tuffs et al. 2004; Mollenhoff et 
al. 2006), but in the case of strong infrared-selected lenses 
with magnifications /i > 10, differential magnification is so 
strong as to make integrated SEDs completely unreliable in 
estimating relative bolometric fractions in lensed galaxies, 
unless a lens model is available to determine the differen- 
tial magnification (c.f. e.g. Granato, Danese & Franceschini 
1996, Rowan-Robinson 2000, Farrah et al. 2002, Verma et 
al. 2002, Efstathiou 2006). Nevertheless, for more moderate 
lensed systems (2 < \i < 5), the underlying bolometric frac- 
tions can still be reliably determined even when neglecting 
the effects of differential magnification. 

These results are robust to the assumed configuration of 
the components in the background source. The background 
source would have to be broadly homogeneous in order to 
avoid strong differential magnification effects. One way to 
do this would be to argue that the number of GMC regions 
is much larger, perhaps numbering several tens to resemble 
those in the Milky Way or M31 (e.g. Sheth et al. 2008). 
However, there is little evidence that star formation in high- 
redshift infrared-luminous galaxies so closely resembles z = 
spiral discs. Furthermore, an active nucleus will inevitably 
create a mid/FIR colour gradient in the background source. 

Despite the lens magnification, the prospects are excel- 
lent for using e.g. the [Cil]/FIR vs. CO/FIR line diagnostic 
diagram (e.g. Stacey et al. 2010, Valtchanov et al. 2011) 
for determining the physical conditions in the emission line 
gas. The line ratios of HCN/HNC and [Oii]/[Cii] appear to 
be particularly insensitive to differential magnification. The 
differential magnification effects on the H2O lines can be ne- 
glected for moderate magnifications (2 < fi < 5) but not for 
higher magnifications. Other line diagnostics are nonetheless 
prone to differential magnification, such as the CO ladder, 
even for low magnification systems (/x > 2). For example, the 
line ratio CO( J = 6 — 5)/CO( J =1 — 0) has a la dispersion 
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Figure 12. The monochromatic magnification factors compared 
to the magnifications of specific components, for an example 
source at z = 2 and a lens at z = 1. Note that the AGN is 
well-traced by rest-frame mid-infrared photometry, and cirrus by 
photometry on the Rayleigh- Jeans tail, but that there is no sin- 
gle monochromatic tracer for star forming regions. Measuring this 
component, and decoupling it from the extended cirrus compo- 
nent, requires resolved imaging in the submillimetre or millimetre 
wavelengths. 



L>L 



600 
500 
400 
300 
200 
100 


10000 
8000 
6000 
4000 
2000 





2<fi<5 





0.2 0.4 0.6 0.8 
U\gn/ L b0 | 



Figure 13. Histograms of simulated AGN bolometric fractions 
for z = 4 lensed galaxies, selected at an observed wavelength of 
500 /im from a 100 /im rest-frame Schechter function with a faint- 
end slope of —1.46. Samples are shown brigher than L* (left) 
or 10L* (right). The underlying bolometric fraction is 0.3. The 
moderate magnification (2 < \i < 5) and strong magnification 
(/i > 10) cases are shown. Qualitatively similar results are found 
at other source redshifts, and the results are only weakly depen- 
dent on the faint-end slope of the luminosity function. 
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of rsj 30% due purely to the effects of differential magnifi- 
cation for all the simulations presented in this paper. One 
might try to avoid these strong differential magnification ef- 
fects by arguing the CO transitions are spatially co- located, 
but this is unlikely to be the case. Even in local star- forming 
galaxies, the low J transitions of CO are more spatially ex- 
tended than higher transitions (e.g. Paglione et al. 2004, 
Sakamoto et al. 2011), and evidence is accumulating of sim- 
ilar trends in high-redshift infrared-luminous galaxies (e.g. 
Ivison et al. 2011). 

There is therefore a set of photometric and spectro- 
scopic diagnostics that do not require a detailed lens mag- 
nification model, and others for which such a model is ben- 
eficial or even essential. Snapshot observations with the At- 
acama Large Millimetre Array (ALMA), e-Merlin, the Ex- 
panded Very Large Array (eVLA), the Plateau de Bure In- 
terferometer (PdBI) and/or the Combined Array for Re- 
search in Millimeter- Wave Astronomy (CARMA) of lens 
candidates from Herschel, Planck, SPT or SCUBA-2 would 
easily establish whether a system is in a low-magnification 
(e.g. \i c± 2) or high-magnification (e.g. \i > 10) regime, 
because the image separation would be enough to con- 
strain the lens magnification. Ground-based optical or near- 
infrared imaging may also assist, but the very red colours 
of some lensed submillimetre galaxies (e.g. Hopwood et al. 
2011) suggest that HST near- infrared observations will of- 
ten be required to detect the background source. High- 
resolution submillimetre/millimetre-wave interferometry is 
also the only way to decouple the lensed star-forming knots 
from the lensed extended submillimetre emission (Figs. [11] 
and !12|) . giving an additional clear rationale for ALMA and 
SMA follow-ups. 

This paper has focussed on FIR to millimetre-wave 
emission line diagnostics, but similar considerations may 
be applicable to rest-frame optical and ultraviolet lines. 
The presence of Lyman a emission in some heavily-obscured 
submillimetre-selected galaxies (e.g. Chapman et al. 2005) 
necessarily requires highly anisotropic obscuration, perhaps 
through the evacuation of cavities via supernova- driven 
winds. Differential magnification may clearly affect the oc- 
curence of nebular emission lines in lensed infrared-luminous 
galaxies. If the nebular lines are close to or coincident with 
the GMC components in the galaxy, then there may be 
examples of submillimetre/millimetre-wave-selected lensed 
infrared-luminous galaxies in which the UV/optical emis- 
sion lines are more easily detectable. A similar considera- 
tion may explain the strong narrow AGN emission lines in 
IRAS FSC 10214+4724, and differential magnification may 
also help explain their unusual line ratios and widths (e.g. 
Serjeant et al. 1998). However, it does not necessarily fol- 
low that close proximity leads to enhancement. The model 
GMCs have a 50 pc radius, embedded in a warm CO knot 
of 400 pc radius, yet the high- J CO lines do not necessar- 
ily share the total magnification boost of the GMCs. Simi- 
larly, extended Lyman a nebulosity may not necessarily be 
as boosted as the continuum sources at the selection wave- 
length. 

If the images of the background source are resolved and 
a lens mass model can be constructed, one might alterna- 
tively choose to regard differential magnification as an asset, 
rather than a problem. Lensing conserves surface brightness, 
so high magnification regions also have high angular mag- 



nification. The presence of caustics in the foreground lens 
magnification could be used to make uniquely high angular 
resolution reconstructions of active nuclei or star-forming 
regions in infrared-luminous galaxies. 



5 CONCLUSIONS 

In infrared-selected gravitational lenses of any magnifica- 
tion fi > 2, the CO ladder is strongly distorted by differ- 
ential magnification. The bolometric fractions of the FIR 
[Cll] and [Ol] lines are also affected by differential magnifi- 
cation, but the distortion is not typically sufficient to alter 
the physical interpretation. Similarly, differential magnifica- 
tion should not typically move a galaxy more than a factor 
of about two from its position on the proposed galaxy "main 
sequence". For moderate-magnification lenses (2 < \i < 5), 
the FIR H 2 0, HCN and HNC lines are broadly unaffected 
by differential magnification, and the broadband SED-based 
decompositions of bolometric contributions should be reli- 
able. However, for strong lenses (/j, > 10) in the absence of a 
foreground mass model, the H2O lines are significantly dis- 
torted, and unless a detailed lens model can be constructed, 
the apparent bolometric contributions of AGN, GMCs and 
cirrus derived from the SEDs are so unreliable as to be use- 
less. 



6 NOTE ADDED IN PROOF 

After this paper was accepted, Hezaweh et al. (2012) 
presented a complementary study of differential magni- 
fication for a source population with varying physical 
sizes. They found that highly- magnified populations over- 
represent compact physical components in good agreement 
with Section 3.4 of this paper. 
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